Preprint typeset in JHEP style - PAPER VERSION 



IPPP/09/70 
DCPT/09/140 
Cavendish-HEP-09/17 



Central Exclusive \, Meson Production at the 
Tevatron Revisited 



L.A. Harland-Lang 1 , V.A. Khoze 2 3 , M.G. Ryskin 2 4 , W.J. Stirling 



1.2 



1 Cavendish Laboratory, University of Cambridge, J.J. Thomson Avenue, Cambridge, 
CBS ORE, UK 

2 Department of Physics and Institute for Particle Physics Phenomenology, University 
of Durham, DEI 3LE, UK 

3 School of Physics & Astronomy, University of Manchester, Manchester MIS 9PL, UK 

4 Petersburg Nuclear Physics Institute, Gatchina, St. Petersburg, 188300, Russia 



Abstract: Motivated by the recent experimental observation of exclusive Xc events at 
the Tevatron, we revisit earlier studies of central exclusive scalar XcO meson production, 
before generalising the existing formalism to include Xci and Xci mesons. Although XcO 
production was previously assumed to be dominant, we find that the Xci an d Xc2 rates for 
the experimentally considered Xc — ► J/i^l — > M + A*~7 decay process are in fact comparable 
to the XcO rate. We have developed a new Monte Carlo event generator, SuperCHIC, 
which models the central exclusive production of the three Xc states via this decay chain, 
and have explored possible ways of distinguishing them, given that their mass differences 
are not resolvable within the current experimental set-up. Although we find that the 
severity of current experimental cuts appears to preclude this, the acceptance does not 
change crucially between the three states and so our conclusions regarding the overall rates 
remain unchanged. This therefore raises the interesting possibility that exclusive Xci and 
Xc2 production has already been observed at the Tevatron. 



1. Introduction 



The measurement of central exclusive production (CEP) processes in high-energy proton 
- (anti) proton collisions represents a very promising way to study the properties of new 
particles, from exotic hadrons to the Higgs boson, see for example Refs. [1] - [7]. 
The CEP of an object A may be written in the form 

pp(p) -► p + A + p(p), 

where + signs are used to denote the presence of large rapidity gaps. An attractive ad- 
vantage of these reactions is that they provide an especially clean environment in which 
to measure the nature and quantum numbers (in particular, the spin and parity) of new 
states, see for example Refs. [4, 8, 9]. A topical example is the CEP of the Higgs bo- 
son [10] - [14]. This provides a novel and promising way to study in detail the Higgs sector 
at the LHC and gives a strong motivation for the addition of near-beam proton detectors 
to enhance the discovery and physics potential of the ATLAS and CMS detectors at the 
LHC [15] - [17]. 

Recently, exclusive diffractive processes pp — > p + A+p have been successfully observed 
by CDF Collaboration at the Tevatron, where A = 77 [18], dijet [19] or \c [20]. 1 As the 
sketch in Fig. 1(a) indicates, these processes are driven by the same mechanism as exclusive 
Higgs (or other new object) production at the LHC, but have much larger cross sections. 
They can therefore serve as "standard candles", see [12, 23], which allow us to check the 
theoretical predictions for the CEP of new physics signals by measurements made at the 
Tevatron. Moreover, the observed rates of all three CEP processes measured by the CDF 
collaboration are in broad agreement with theoretical expectations [3, 5, 12, 23], which 
lends credence to the predictions for exclusive Higgs production at the LHC. 

Among the CEP processes measured at the Tevatron, the double-diffractive production 
of C-even, heavy quarkonia (% c ) states plays a special role [5] (see also [24] - [29]). First, 
as is well known, heavy quarkonium production provides a valuable tool to test the ideas 
and methods of the QCD physics of bound states, such as effective field theories, lattice 
QCD, NRQCD, etc. (see, for example, Refs. [30, 31] for theoretical reviews). Second, 
Xc production exhibits characteristic features, based on Regge theory, that depend on the 
particle spin and parity J p . We discuss these in Section 2 below. 

A potential problem with Xc production as a "standard candle" for Higgs produc- 
tion is that it is far from clear that the purely perturbative approach of Refs. [10, 12] 
(as exemplified by Fig. 1(a)) is valid. In particular, the estimates in Ref. [5] assume a 
perturbative contribution coming from integrating round the gluon loop in Fig. 1(a) for 
Q± > 0.85 GeV. Due to such low scales, strictly speaking one cannot guarantee that the 
accuracy in the perturbative predictions made in this way for Xc CEP is better than a 
factor of 4-5 up or down. In addition, there may be important non-perturbative contri- 
butions, traditionally modelled by the Pomeron-Pomeron process shown in Fig. 1(b), but 

1 For a recent review see [21]. More CDF exclusive data on 77 and \c production may be available in 
the near future [22]. 
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Figure 1: Schematic diagrams for CEP of a system A within the approach of Refs. [10, 12] and 
[32] - [35]. The integration over the loop momentum Q± in diagram (a) results in a J z = selection 
rule [3], where J z is the projection of the total angular momentum along the proton beam axis. It 
is also necessary to compute the probability, S 2 , that the rapidity gaps survive soft ([32, 33]) and 
semi-hard ([34] - [36]) rescattering; these two possible types of unitarity (or absorptive) corrections 
are exemplified in diagrams (c) and (d) respectively, where the dashed lines represent Pomeron 
exchanges (as in version (b) of diagram (a)). 



here again there is significant model dependence and so any predictions for this component 
also come with large uncertainties. We choose to take a pragmatic approach, in which we 
base our analysis on the perturbative contribution only, but at the same time we consider 
which features of the perturbative contribution (for example, the relative contributions of 
the various J p states, distributions of final state particles, etc.) are likely to be shared 
by the non-perturbative contribution. Independent of the exact details of the production 
mechanism, the CEP of x states provides a valuable check on the important ingredients of 
the physics of Pomeron-Pomeron fusion. 

Another important issue related to the CEP of \c states at the Tevatron is that this 
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process allows us to test the role of the so-called enhanced absorptive corrections (see 
Fig. 1(d)), which break soft-hard factorization, see for example [35] - [36]. As shown in 
[14, 35], there is a hierarchy in the value of the rapidity gap survival factor S due to 
enhanced absorption, S cn h, 

S™ C (M H > 100 GeV) > ^r ron (77;^?>5GeV)> Ch atr ° n (Xc), 

which reflects the size of the various rapidity gaps (s/M 2 ) of the different exclusive pro- 
cesses. The very fact that 77 and Xc events have been observed at the Tevatron in rea- 
sonable agreement with theoretical expectations 2 confirms that there is no danger that 
enhanced absorption will strongly reduce the exclusive SM Higgs signal at the LHC. 

When interpreting the results of [20] in terms of the production of a particular \c 
state there is one important point to bear in mind. While the P and C parities are 
unambiguously defined by the fusion mechanism (see, for example, [12]) the spin assignment 
requires special care. As discussed in Refs. [3, 5], central exclusive Xc production should 
be dominated by the Xco(0 ++ ) state. This is because Xci(l ++ ) and Xc2(2 ++ ) production 
is strongly suppressed: the former due to the Landau- Yang theorem [38] for on-mass-shell 
gluons and the latter because in the non-relativistic approximation the Xc2(2 ++ ) meson 
cannot be produced in the J z = state, which dominates CEP for forward outgoing 
protons [3]. 

Recall, however, that the experimental observation [20] of exclusive Xc production is 
based on the decay chain Xc — * J 1^1 ~^ ^ + ^1- The observed 65 ± 10 signal events have a 
limited M(J/ipj) resolution and are collected in a restricted area of final state kinematics 
(due to cuts and event selection criteria) . In order to determine the Xc yield the dominance 
of XcO production is assumed , and the CHIC Monte Carlo 3 , based on the XcO — ► J/V' + 7 
decay, is used for conversion of the observed events into the cross section. However at 
the present time we cannot rule out the possibility that, under the conditions of the CDF 
experiment, higher spin Xc states (1 ++ , 2 ++ ) contribute to the observed J/ip + 7 signal. 4 
As is correctly pointed out in Ref. [29], the strong suppression of 1 ++ central production 
can be compensated by its much higher branching fraction to the J /ip + 7 final state. 
We show below (see also [14]) that this can also be true for Xc2(2 ++ ) CEP. Explicitly, 
the XcO: Xci, Xc2 branching fractions to J/ipj are 0.011, 0.34 and 0.19 respectively [37]. 
There is another factor leading to a further rebalance between the relative contributions of 
different Xc spin states. As discussed in [4, 14], the eikonal survival factor, S^, is larger 
f° r Xci and Xc2 since, due to their spin structure, they are produced more peripherally. 5 

The simultaneous presence of several Xc states clearly requires a more comprehensive 
analysis, including a new Monte Carlo programme, allowing for production and decay of 

2 In the Xc case the agreement becomes especially striking after taking into account the revised value of 
the total XcO width which has been reduced by a factor 1.4 [37] as compared to the value in the Review of 
Particle Properties (2002) used in [5]. 

3 CHIC is a publicly available Monte Carlo implementation of the \c0 analysis of Ref. [5]. 
4 This applies not only to % C (1P) states, but also to possible higher excitations Xc(nP). 

5 *S en h is largely independent of the \ s P m assignment. Note also that the relative number of events 
where the forward protons dissociate is larger for Xci and \c2 than for XcO- 
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the higher-spin states. This is the main topic of the current paper, i.e. we extend the 
analysis of Refs. [3, 5] to include the detailed study of Xd and Xc2 exclusive production. 
Special attention is paid to the role of absorptive corrections, which significantly affect the 
predicted rates, see [4, 5, 14]. 

The paper is organized as follows. In Section 2 we review the general expectations 
for the J p properties of Xc production based on Regge theory. In Section 3 we discuss in 
detail the perturbative approach to the calculation of Xc CEP, paying particular attention 
to the differences between XcO, Xd and Xc2 production and to the uncertainties in the 
predictions. We have implemented our calculations in a new Monte Carlo event generator 
- SuperCHIC - which is described in Section 4. In Section 5 we present numerical results for 
XcO, Xd and Xc2 CEP at the Tevatron. We discuss the impact of the survival factors on the 
CEP cross sections, and comment on the size of possible non-perturbative contributions. 
We investigate to what extent kinematical distributions of the final-state particles can be 
used to distinguish the three Xc spin states. In Section 6 we summarise our conclusions and 
comment on possible future developments in the study of the CEP of quarkonia states at the 
Tevatron and LHC. Some additional calculational details are presented in two Appendices. 



2. General expectations from Regge theory 

As discussed in [4], the central diffractive production of meson states (see Fig. 1(b)) has 
characteristic features that depend on the particle spin and parity J p , which follow from the 
general principles of Regge theory. Let us first recall particular examples of bare Pomeron- 
Pomeron vertices for the spin-parity J p of particle h in the case of low transverse momenta 
Pi,2 ± °f the outgoing protons. 

(a) J p (h) = 0+ 

For a scalar particle h, the vertex coupling is simply 6 

9pp = /o+(Pi X 'PL'Pu-P2±). 

where / + is a function of the displayed scalar variables. When pf ± or p| ± — ► 0, this 
function, in general, tends to some constant f s . Further information on the structure 
of this function requires extra dynamical input. In particular, within the perturbative 
framework (for Q 2 ± » Pi :2± ), fo+ is almost independent of Pi,2 ± , and hence the bare 
cross section, a + , is essentially independent of the azimuthal angle <j> between the 
outgoing protons 

d(T r 0+ , oc constant (0), (2.2) 

where — — P?2 ± - We note, however, that for Xc CEP the low Q± scale we are 
considering means that the inequality Q x Pi 2 x is not completely valid, and we 
therefore expect some deviation from the constant behaviour of (2.2). 



°We use boldface type to denote spatial three-vectors. 
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(b) J p (h) = o- 

For the central production of a pseudoscalar particle, the bare vertex factor takes the 
form 

9pp = /o-(Pi ± 'P2 ± 'Pi±-P2±) (Pi± x P2 X ) • n, (2.3) 

where n is the unit vector in the direction of the colliding hadrons (in the c.m.s.). 
Due to the identity of the Pomerons, the function / - should be symmetric under 
the interchange 1 <-> 2. It follows from (2.3) that in the pseudoscalar case the bare 
cross section, cr -, should behave for small |i 1,2 1 as 

W K 1(111(21 si,lV (24) 

An immediate consequence of (2.4) is that pseudoscalar production is forbidden when 
the protons scatter at zero angle. 

(c) J p (h) = 1+ 

For the production of an axial vector state the bare Pomeron-Pomeron fusion vertex 
factor can be written as 

n A „ (ti -t 2 )([ Pu x P2 J -e) [K x n] • e 

9 PP ~a A=0 — 2 + oa=i (2.5) 

with 

K = pi- P2 . (2.6) 

Here M and e are the mass and polarization vector of the centrally produced 1 ++ 
state, and the vertex functions a^=o,i correspond to axial meson production with 
helicities A = 0, 1 in the target rest frame (where the longitudinal component of the 
axial meson momentum is much larger than its transverse component). Analogously 
to / + and f - in the previous cases, the functions a A= o,i may depend on pf ± , p 2± 
and (pi ± • P2 ± ), and are symmetric under the 1 <-> 2 interchange. 

It follows from (2.5) that the bare amplitude tends to zero at low K±, in particular 
when both protons scatter at zero angle. Another important consequence of (2.5) is 
that at low |tx,2 1 the axial meson should be produced dominantly in the helicity-one 
state. As already mentioned in [4], the general structure of the axial vertex gp P , 
given by Eq. (2.5), coincides with that found using a non-conserved vector current 
model [8], which gives a good description of the experimental data on ^(1285) and 
/i(1420) CEP by the WA102 Collaboration [39]; for a review see Ref. [7]. 

(d) J p (h) = 2+ 

For a tensor particle h, the bare vertex function gj, p is not constrained by Regge 
theory alone. However, as we have already mentioned, within the perturbative ap- 
proach of Refs. [5, 12] the forward CEP of non-relativistic heavy 2 ++ quarkonium 
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should be strongly reduced because of the suppression of the 2 ++ — > 2g transition for 
the J z = on-mass-shell two-gluon state. 7 




Figure 2: Impact of the absorptive corrections on the distribution (in arbitrary units) of the 
difference in azimuthal angle of the outgoing protons for the CEP of various J p y c states at the 
LHC, using the two channel eikonal model of Ref. [32]. The solid (dashed) lines are the distributions 
including (excluding) the survival factor. For completeness we also show the result for pseudoscalar 
rj c production. 



Finally, we recall that, as discussed in [4], the absorptive corrections arising from the 
multi-Pomeron exchanges modify the distributions over (j). This is because the absorption 
depends on the distribution in impact parameter b space, which in turn leads to a char- 
acteristic dependence of the survival factor (mainly S^ ik ) on the azimuthal angle between 
the outgoing protons. This effect was discussed in detail in [40]. 

To demonstrate how the absorptive corrections may affect the angular distributions 
between the outgoing protons we have used a simple two channel eikonal model [32]. We 
show in Fig. 2 the results for ^fs = 14 TeV. An analogous dependence of Sr ik on the 
transverse momentum of the centrally produced meson is shown for the Tevatron energy 
y/s = 1.96 TeV in Fig. 3. 

7 We reconfirm the conclusion of [3] that the relativistic corrections to the Xc,f>(2 ++ ) — > 2g transition are 
numerically small. 
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Figure 3: Dependence of the survival factor S^ ik on the transverse momentum of centrally produced 
Xc mesons at the Tevatron, using the two channel eikonal model of Ref . [32] . Also shown (dashed 
line) is the survival factor corresponding to r\ c production. 



3. Central Exclusive Xc production: perturbative framework 

To calculate the perturbative contribution to the central exclusive \c production process 
we use the formalism of Refs. [3, 4, 10]. The amplitude is described by the diagram shown 
in Fig. 1(a), where the hard subprocess gg — > Xc is initiated by gluon-gluon fusion and the 
second t-channel gluon is needed to screen the colour flow across the rapidity gap intervals. 
We can write the Born amplitude in the factorised form [5, 6] (see Fig. 4): 

T = tt 2 [ —2— d ^tji. 72 • fg{xi,x' 1 ,Ql,ti 2 ;t 1 )fg(x2,x' 2 ,Ql,n 2 ;t 2 ) , (3.1) 

where Vj is the colour- averaged, normalised sub-amplitude for the gg — >■ Xcj process: 

v ' = ^tf^r[T,**rtrt v $- (3.2) 

C a,b 

Here a and b are colour indices and Nc = 3. The amplitude represents the coupling of 
two gluons to the Xc state being considered: the procedure for calculating this is outlined 
below. The / s 's in (3.1) are the skewed unintegrated gluon densities of the proton at the 
hard scale fx, taken typically to be of the order of the produced massive state, i.e. M x /2 in 
this case, and only one transverse momentum scale is taken into account by the prescription 

Qi = min{Q_L, |(Q ± - pi ± )|} 

Q 2 = min{Q ± ,|(Q ± + P2± )|}. (3.3) 
The longitudinal momentum fractions carried by the gluons satisfy 
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Figure 4: The perturbative mechanism for the exclusive process pp — > p + x + with the eikonal 
and enhanced survival factors shown symbolically. 



Q± 



<C x 



M 



(3.4) 



The t dependence of the / 5 's is not well known, but in the limit that the protons scatter 
at small angles, we can assume a factorization of the form 



f g (x, x', Q 2 ,^ 2 ;t) = f g (x, x', Q 2 , n 2 ) F N (t) , 



(3.5) 



where the i-dependence is isolated in a proton form factor, which we take to have the 
phenomenological form F^{t) = exp(6t/2). In such a case a fit to soft hadronic data [32] 
gives b ~ 4GeV -2 , which is also consistent with extracting Fj^(t) from the t-dependence 
in 'elastic' J/ip photoproduction [41]. We will therefore use b = 4GeV~ 2 as our value for 
the slope parameter throughout. 

In the kinematic region specified by (3.4), the skewed unintegrated densities are given 
in terms of the conventional (integrated) densities g(x,Q 2 ). To single log accuracy, we 
have 8 

f g (x, x', Qlv 2 ) = R 9q^q2[M*, Q 2 )^T g (Qlv 2 )] , (3.6) 

where T g is the usual Sudakov survival factor which ensures that the active gluon does 
not emit additional real partons in the course of the evolution up to the hard scale ti, so 
that the rapidity gaps survive. R g is the ratio of the skewed x' <C x unintegrated gluon 
distribution to the conventional diagonal density g(x,Q 2 ). For x <C 1 it is completely 
determined [43] . The explicit form for T g is given by resumming the virtual contributions 
to the DGLAP equation. It is given by 



T g (Q 



2 ^ 2 ) 



exp^- 



dk\ a s (k 2 ± ) 



2tt 



l-A 



dz 



(3.7) 



Here, as in [5], we go beyond the collinear approximation and in the T factor we resum not 
just the single collinear logarithms, but the single soft ln(l — z) terms as well. To a good 
approximation, this can be achieved by taking the upper limit of the z integration in (3.7) 



3 In actual calculations, we use a more precise phenomenological form given by Eq. (26) of [42]. 
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to be 



A = 



(3.8) 



k L + 0.62M X 



Returning to the gg — > \c amplitude, we note that the original extension of the CEP 
formalism to XcO production in [5] was achieved in direct analogy to the Higgs case, that is 
by assuming that the XcO coupled to the gluons as a pure scalar with any effects from its 
internal structure neglected. We now wish to go beyond this approximation and model the 
internal structure of the \ c meson for all three J states and in particular their coupling to 
two gluons. This is done by a simple extension of the calculation of [44], where the coupling 
of 3 Pj quarkonium states to two off-mass-shell photons is considered: as the gluons are in 
a colour singlet state the only difference will be constant prefactors resulting from colour 
algebra. We will simply state the results for the three amplitudes, leaving the derivation 
to Appendix A: 



where qi ± = Qj_ — p± ± and q2 ± = —Q± — P2 ± - The amplitudes are normalised as in (3.2) 
and the qi ± are 4-vectors with qf ± = — q? < throughout. 9 Considering first the XcO 
vertex, in the <C M 2 limit (which is true to an acceptable degree of accuracy) we 
expect the internal structure of the XcO to be unimportant, and therefore to recover the 
previous result of [5]. We find 




V 



(3.11) 



(3.10) 



(3.9) 



Vq 



487ra 5 ^(0) 




(3.12) 




Making use of the standard NRQCD result (see for example Refs. [45, 46, 47]), 



r(x c0 -<75) = 96^||^(0)| 2 , 



(3.13) 



9 Four- vector scalar products are denoted by (pq). 
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we find 10 

' °' = M3 ^U^J . ( 3 - 14 ) 

which has the same form and normalisation as the previous result, as it must do. We will 
take this large M x limit throughout. Turning now to the Xd vertex, we can immediately see 
that it vanishes for on-shell gluons, that is when qf = q? ± = 0, as dictated by the Landau- 
Yang theorem (see Section 2). Furthermore, in the forward limit we have qi± = —q2± = Q± 
and so 

r 3 
— ( 

2 



V ^-\/^cM x Ql, (3.15) 



Aic 



Vi - —Q 2 ± Pi,uP2, a Q±^ef , (3.16) 
y 2 _y^( s Q ±M Q ±a + 2Qlp lllP2a )e*/ a . (3.17) 

We see that V\ is odd in Q±, and will therefore vanish upon the loop integration (3.1) over 
Q_i_. For V2 we make use of the identity 

j d 2 Q ± Q ±ll Q ±a = I J dQlQlg^ , (3.18) 

where g^ a , the transverse part of the metric, can be written in the covariant form 

2 

9la = 9fia ~ ~iPllJ.P2a +PlaP2n) ■ (3.19) 

We then find V 2 oc which vanishes due to the tracelessness of the \2 polarization 
tensor (A. 18). We see that, as expected, the Xc2 and Xd production amplitudes vanish in 
the forward limit, and we will therefore expect the corresponding rates to be suppressed 
relative to XcO production, via the integration over the proton form factor e bti ~ e~ bp i->- 
which suppresses large pf ± values. In fact we can give a very rough estimate for the level 
of suppression we will expect. Squaring and summing over polarization states gives 

IT/ I 2 . IT/ I 2 . IT/ |2 1 . (P±) . (Pi-) (o orA 

|Vo| • |Vi| . \V 2 \ ~ 1 . ■ — , 2 • (3.20) 



X 



Note that this result, as well as the amplitudes of (3.9 - 3.11), is also applicable to the CEP 



Analogously to Ref. [5], we assume the same NLO correction for the gg — > x vertex as for the X ^ 99 
width, which can be valid only within a certain approximation. Moreover, as has been known for some 
time in the P— wave case (see for example Ref. [48]), the NNLO and higher-order radiative corrections 
to the X ~> 99 transition are expected to be numerically quite large, and this would result in further 
uncertainties in the theoretical expectations. Recall that these corrections are not universal and depend on 
the spin-parity assignment of the P— wave states. In particular, it is known that in the xo case the dominant 
part of the NLO correction comes from the {in) 2 term originating from the Sudakov-like double logarithm 
a s ln 2 (q 2 /M 2 ), when the imaginary part of the logarithm (= — in) is squared. This double logarithm, 
and correspondingly the {in)' 2 contributions, are absent for the case of the xij where the amplitude for 
on-mass-shell {q 2 — 0) gluons vanishes due to the Landau- Yang theorem. 
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of Xb mesons. The factor of <p^_) comes from integrating over the assumed exponential 
form of the proton vertex, 

<pl> = / dple- b ^ = I = \ GeV 2 . (3.21) 

If for simplicity we assume (Qj_) ~ 1.5 GeV 2 and 

M 2 » 10 GeV 2 , we obtain 11 

|Vb| 2 :|Vi| 2 :|V2| 2 ~l:^:^. (3.22) 

While it is clear that we will have a quite sizeable suppression of the Xci an d Xc2 CEP cross 
sections, these values are of course only very rough estimates, and an explicit calculation 
is required to confirm them. 

We stress again that the legitimacy of the extension of the purely perturbative QCD 
treatment for central exclusive Higgs production to the Xc case is somewhat questionable. 
For Higgs production the hard scale fi is set by M^/2, and so we expect that a reliable 
calculation within perturbative QCD can be performed. In particular, the Sudakov factor 
leads to an IR stable result, with only a small contribution to the cross section from the 
region of Q± below ~ 1 GeV (although this is not to say the calculation does not come 
with significant uncertainties). However, in the case of Xc production, where the 'hard' 
scale is ~ 1 GeV, we expect and find that a significant part of the cross section comes 
from the IR unstable low Q± region. It might seem then, that despite the attractions of 
considering central exclusive Xc production, any attempt to calculate a reliable cross section 
within the perturbative QCD framework is unlikely to succeed. This is not the case: the 
philosophy we take is that, even if the purely perturbative calculation is not IR stable, we 
should expect a smooth matching between the perturbative regime and the 'soft' regime 
to which we can apply a non-perturbative Regge model. This was done in [5] for the XcO, 
with a simple model for the Pomeron invoked, and the non-perturbative and perturbative 
contributions were found to be of a similar size, which gives justification for the inclusion 
of a perturbative contribution to Xc CEP. On the other hand, there is much uncertainty 
surrounding which non-perturbative models for the Pomeron are most appropriate in this 
context, and so results such as this can only be used as a guide. 

Returning to the perturbative calculation, several comments are in order. As men- 
tioned above, in the case of the Xc we expect a significant proportion of the cross section 
to come from the low Q± region where perturbation theory is not valid, and there will 
correspondingly be a large degree of uncertainty in its predicted value. Given our lack of 
detailed understanding of low Q± non-perturbative gluon dynamics, the best we can do is 
to introduce an infrared cut-off to the Q± loop integral such that we are only considering 

11 In fact, accounting for the values of pi±,p2± in the denominator of (3.1) we obtain a slightly larger 
value of the effective slope 6 eff = b + 0(l/Q 2 ± ) > 4 GeV -2 . Therefore we expect a slightly larger suppression 
of the higher spin, xi and \2, states than that given by (3.22). Note also that for a heavier meson (xb) 
CEP the slope b cB will be smaller than that for the case of \c due to a typically larger values of Q± in the 
integral (3.1). 
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the regime where perturbation theory will be reliable. This can be loosely justified on 
the grounds that our current understanding of these non-perturbative dynamic predicts 
that the low Q± contribution appears to be suppressed [49]. Nevertheless significant un- 
certainties remain, both in the specific choice of cut-off, for which we have rough physical 
guidelines but which inevitably amounts to a subjective decision, and in what contribu- 
tion we would actually expect from the low Q± region, which a cut-off prescription simply 
ignores. 

A further uncertainty that is worth mentioning arises from the skewed PDFs, which we 
can express via (3.6) in terms of the conventional PDFs g(x,Q 2 ), and therefore evaluate. 
Unfortunately, although this is in principal true, there is a large degree of uncertainty in 
the value of the conventional PDFs at the low x and low Q 2 scales we are considering, as 
can be seen in Fig. 5 where four representative PDF sets are plotted at Q 2 = 1.5 GeV 2 . In 
the case of exclusive Xc production at the Tevatron (y/s = 1.96 TeV) we are sampling the 
x ~ 2 x 10 -3 region, for which there is a large uncertainty. Recalling that the final cross 
section depends quartically on the skewed PDFs, this is clearly unacceptable. We can see, 
on the other hand, that as x is increased the uncertainty rapidly decreases and we have 
a reasonable agreement between the sets in the region 0.05 < x < 1. We can therefore 
perform the cross section calculation for a lower c.m.s energy such that the sampled x 
value is in this range, although we must be careful that the corresponding x value is not 
too high, as our initial formula for the skewed PDFs relies on a small x approximation. In 
actual calculations a value of \/s = 60 GeV is chosen, which corresponds to x ~ 0.05. To 
make contact with the experimental c.m.s. energies of the Tevatron we then assume that 
the total cross section exhibits the Regge behaviour 

da OC s ap(.ti)+a P (t2)-2 _ (3^3) 

This gives us a simple way to avoid the large uncertainties of the PDFs in the low x region, 
although our final result will depend on the validity of this Regge assumption and on the 
specific value of the Pomeron intercept ap(0) that we use. 12 

We end this section with a brief review of the uncertainties that are present in our 
calculation. First, we have the uncertainty in our choice of hard scale /j, and the prescription 
(3.3) for the transverse momentum scale Q 2 , with for example the choice [i = Mh/4 in 
the Higgs case giving a quoted enhancement of 30% to the cross section [4]. Further to 
this we have the even more considerable uncertainty coming from the dependence of the 
cross section on the conventional PDFs to the fourth power in the low x and Q 2 region, 
where they are not well determined. We also have the dependence of the final result 
on the non-perturbative, non-universal survival factor, which gives perhaps the largest 
uncertainty in the overall production rate. Differences in the higher-order QCD radiative 
corrections to the gg — ► x vertex could also cause additional uncertainties. Finally we have 
the uncertainty in the low Q± contribution to the perturbative amplitude, which we find 
to be quite large in the Xc case. For more details of these issues we refer the reader to [4]. 



Our procedure is equivalent to assuming a Regge-based extrapolation of the PDFs to small x. 
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Figure 5: LO and NLO PDFs at Q 2 — 1.5 GeV 2 , plotted as function of x. Four representative 
PDF sets are displayed, and a large uncertainty at small x is clear. 

4. SuperCHIC Monte Carlo Generator 

In [5], central exclusive XcO meson production and its subsequent decay to J/tfij — > /i + ^~7 
was modelled using the CHIC Monte Carlo event generator, but we now wish to consider 
the case where higher spin \c states are produced. The CEP of the three Xc(o,l,2) states is 
modelled using a new, more general, Monte Carlo programme, SuperCHIC. 13 This follows 
essentially the same procedure for generating the relevant phase space as the previous 
CHIC MC - optimised to reduce the event weight variation - but with some important 
generalisations included that we will now outline. 

The explicit evaluation of (3.1) is in all cases performed 'offline' from the SuperCHIC 
Monte Carlo event generator: to perform the loop integration for each event would lead to 
an unacceptably large run-time. The skewed PDFs are calculated as outlined in Section 3, 
with the Sudakov factor read in from a grid to minimise run-time, and (3.1) then evaluated 
using standard Monte Carlo techniques. We go beyond the approximation used in [5], 
where the amplitude squared was calculated in the forward limit, with the p± dependence 
isolated in the proton form factors. Such an assumption is not relevant for the Tevatron, 
where the p± of the final state protons is not measured [20], and so we must include 
non- forward effects. This will not only give a more accurate evaluation of the XcO cross 
section, but is also essential in the case of Xd and Xc2 production, where the corresponding 
amplitudes vanish in the forward limit. On the other hand, we still need to perform the 
loop integration separately from the Monte Carlo event generator. Considering first the 
XcO amplitude Aq, this can be achieved by noting that for small p± it must have the form 



13 The extension to, for example, the CEP of pseudo-scalar r/ c or higher excitation \c{nP) states, as well 
as the respective 6-quark states Xb, Vb an d Xb(nP), is planned for future work. 
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that is, there exists a Taylor expansion for Aq formed from all possible scalar combinations 
of the pi ± , the validity of which depends on the suppression in p\ coming from the proton 
form factors. Squaring (4.1) and keeping only the leading terms in p? , we can see that 
this expansion is equivalent to making the replacement (at lowest order in pf ± ) 



e 



M ± _ e -(^-^)p? ± . (4 . 2) 



Thus to a first approximation we expect the inclusion of non-zero p± in the amplitude 
calculation to simply result in a change in the effective slope of the proton form factor. 
This then allows for an easy way to take into account the effect of non-forward protons in 
the amplitude, as we can simply model \Aq\ 2 as a Gaussian, 

e -KP? x+ Pi x )| Ao |2 0(e -C(p? ± +Pi ± ) ) (4 . 3) 

where the slope and the overall normalisation are set by matching the values of (p£ ± ) 
and the integrated cross section, respectively, to those given by the exact expression for 
|^4o| 2 - We note that the exact expression for Aq (and therefore the p x± distribution) will 
depend in general on the azimuthal angle between the outgoing protons, but by choosing to 
model this effect by simple Gaussians in pf ± the resultant azimuthal correlations between 
the protons will not be fully modelled in the Monte Carlo. Thus in the case of the XcO 
any deviation, for example, from the flat behaviour of (2.2) is ignored. On the other hand 
it is clear that we are not currently interested in correctly modelling the p± distributions 
of the outgoing protons, which we recall are not measured at the Tevatron, but only 
those of the centrally produced final state particles, which it is important to know when 
discussing possible methods for distinguishing the three Xc states, and the simple Gaussian 
approximation achieves this to an acceptable degree of accuracy. 14 However, we should be 
careful in our application of this approximation, as no strict Q± 3> p± hierarchy exists for 
the Xc- 15 

The higher J states are more complicated due to their non-trivial Lorentz structure, 
but the basic argument remains the same. We can write the Xd and Xc2 amplitudes, 
omitting the Xc polarization vectors etc. for simplicity, as 

A* oc ( P2± - p 1± re~ bf(pi ± +p ^ )/2 , (4.4) 
A? oc(Kp 1± )^(p 2± )^+2(p 1± p 2± Kp^ e - 6 " (p t+PL)/ 2 , (4.5) 

where we must now square the amplitudes and sum over the relevant Xc polarization states 
before matching the values of the normalisation and slope as before. 

SuperCHIC is a standard MC event generator that calculates the relevant weight for 
each generated event using these effective slopes, which can simply be read in at the 



14 In fact, in the case of the Xc2 it is necessary to include a pi^ • p2 ± term in the fit. We also note 
that, while it is not done here, the complete inclusion of the correct azimuthal correlations, which may be 
relevant for measurements at the LHC with tagged forward protons, remains a possible future extension of 
the Monte Carlo. 

15 Note that such a procedure can be useful in the case of CEP of the Higgs boson [6, 12]. 
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beginning of the run. We give the option of generating events for purely XcO, Xci and Xc2 
production, as well as the (experimentally relevant) option of generating all three states at 
once. It is possible to generate the differential cross section da/dy x at a given rapidity y x 
value or the full cross section over a pre-specified Xc rapidity range, with an approximate 
phenomenological fit invoked for the y x dependence of the cross section. 

Having generated the appropriately weighted central exclusive Xc event, we then gen- 
erate the decay process Xc — > J Wt — * ^ + ^~7 through which exclusive Xc production has 
been observed at the Tevatron. While the isotropic decay of the scalar XcO — ► J VV'T is triv- 
ial, the situation for Xci and Xc2, which have non-trivial polarization states that must be 
accounted for in the relevant decays, is not so simple. The calculation of the different decay 
distributions is, however, relatively straightforward and is outlined in full in Appendix B. 
In all cases, these decays have an integrated weight of unity, while we multiply the overall 
cross section by the relevant branching ratios, taken from [37]. For the Xci, Xc2 and J ftp we 
generate the spin states in the helicity basis: that is, we generate the Xc(J/ip) polarization 
vectors in the Xc(J /VO res t frame, before boosting along the spin quantization (z) axis and 
then rotating in the z — P x (^) plane. Finally, to improve the overall efficiency we allow 
as input the option of specifying values for the experimentally significant cuts on the final 
state pair, in particular the maximum pseudorapidity \r]\ and the minimum p±. All 

kinematic information for the produced particles is calculated and can be read out at the 
end of the run. 



5. Results 

We begin, for the sake of comparison with the previous results of [5], with a calculation of 
the XcO cross section in the forward limit. As in [5], we use GRV94H0 partons [50] through- 
out 16 , and we choose the value of Q± = 0.85 GeV as our infrared cut-off. Combining the 
perturbative and non-perturbative contributions, we find 



= 80nb. (5.1) 

yx=o 



Once we have corrected for the revised PDG value for the total XcO width (which has 
decreased by a factor ~ 1-4), as well as the slightly revised value of the survival factor 
that we use, we find that this is in good agreement with the previous result quoted in [5]. 
Moreover, this is also in excellent agreement with the experimental value from the CDF 
collaboration [20]: 



dcr cxp 



'Xc 



= (76±14)nb. (5.2) 

Vx=0 



We emphasise that this value was assumed, rather than being observed, to correspond to 
XcO production, as the mass difference between the three Xc states was not resolvable within 
the experimental set-up. 



16 As already mentioned, the GRV94HO gluon PDF is very consistent with the more recent MSTW2008 
and CTEQ LO and NLO gluon PDFs in the 0.05 < x < 1 region. 
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We now consider the effect of non-forward protons (p± 7^ 0) on this cross section (see 
also [4, 6]). By directly integrating over the p± dependent amplitude squared and fitting 
the resultant p^ ± distribution according to (4.3), we find 

bf = 6.6GeV" 2 , (5.3) 

that is, a steepening in the effective slope of the proton form factor 17 . The effect of this is 
shown in Fig. 6, where we have plotted the differential cross section as a function of p±(xco) 
in the forward and non-forward limits discussed above. As a check, we plot both the exact 
and our fitted results for the non- forward limit: we can see that the match is sufficiently 
accurate for our purposes. As expected, the steeper slope corresponds to p±(xco) being 
more sharply peaked at low p± values. This effect is enhanced by the prescription (3.3) 
for the argument Q\ of the skewed PDFs. In particular, we have introduced an infrared 
cut-off to avoid a contribution from the low Q± domain, where perturbation theory is not 
valid. However we can see that (3.3) will have the average effect of pushing Qi into this 
'soft' region which we are neglecting, and therefore lower values of p± will be favoured, 
leading to a steepening of the effective slope. On the other hand, we could, instead of 
taking the minimum of the two gluon transverse momenta, take their average, which gives 
&Q ff « 5.6 GeV -2 . It is clear then that we expect some steepening in the effective slope, 
but the exact amount is very much tied up in the overall uncertainties of the calculation. 
Naively, we might expect this to lead to a factor of ~3 decrease in the production cross 
section, via the integration over the transverse momenta pi ± , however we recall (see, for 
example, [5, 14, 35]) that the total cross section depends on the ratio S 2 /b 2 , which depends 
only weakly on b 2 ; that is the reduction in the cross section caused by the increased 
slope is largely compensated by an increase in the survival factor for the more peripheral 
interaction, leading to only a small overall decrease in the XcO rate. We next consider the 
Xd and Xc2 cross sections, which are calculated following the procedure outlined in Section 
4. Directly integrating over the p± dependent amplitude squared and fitting the resultant 
p 2 x distribution according to (4.4) and (4.5) for Xd and Xc2 CEP, respectively, we find the 
following values for the effective slopes: 

6f = 4.6GeV- 2 , (5.4) 
bf = 5.9GeV" 2 , (5.5) 

We can see from Fig. 3 that while the XcO and Xd eikonal survival factors are approximately 
constant, the Xc2 eikonal survival factor has a strong dependence on the p± of the Xc2, which 
we read in from a grid in the Monte Carlo. Making use of (5.3), (5.4) and (5.5) we find 

Seik(Xo) : S 2 ik ( X i) : (5 2 k ( X2 )) ~ 0.098 : 0.15 : 0.22 . (5.6) 

We emphasise that the gap survival factor S 2 depends on the effective slope b eS of the bare 

17 Recall that here & eff /2 is the slope of the 'bare' amplitude before screening effects are included. The 
inclusion of absorptive effects (i.e. the gap survival factor) will further enlarge the px-slope of the experi- 
mentally observed forward proton distributions. 
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Figure 6: Differential cross section (in arbitrary units) as a function of Xco P± as a result of exactly 
calculating the p± dependent Xco amplitude ('exact'), setting p± — in the initial amplitude 
calculation ('forward approximation') and approximating the non-forward effects by an effective 
slope parameter b e g ('non-forward approximation'). 



(non-screened) CEP amplitude which, in turn, depends explicitly on the p± of the outgoing 
protons, see [6]. In particular, at b = 4 GeV -2 we obtain S 2 ik = 0.046, while in the case 
of CEP of the XcOi where b® s = 6.6 GeV -2 , we obtain the value 5 2 ik = 0.098, i.e. a factor 
of two larger. As discussed above, such an increase in S 2 largely compensates the decrease 
in the CEP cross section caused by a smaller phase space in p± occupied by the final state 
protons (due to a larger 6 eff ). Recalling that the CEP event rate depends on the ratio 
S 2 /b 2 (rather than on S* 2 ) [5], this can serve as a warning regarding the conclusions made 
in Refs. [51, 52] about the CEP rates based solely on the evaluations of S 2 for protons with 
p± = 0. In addition, we have already shown in Section 3 that neglecting the p± / effects 
in the structure of the hard production subprocess can lead to a significant underestimate 
in the production rate of states with J p other than + . Returning to the calculation of the 
higher spin Xc cross sections we find, by including the relevant branching ratios (evaluated 
at rapidity y x = in all cases), 

"T^O , pert rXl , pert rX2 pe rt 

4^^- ■ -4^^ ■ /* "1 ■ 0-8 : 0.6 . (5.7) 

Ttot dy x r*\ dy x T% t dy x 

Thus, within the perturbative framework, the expected contributions of the three Xc states 
to the Tevatron data are of comparable size, despite the initial suppression in the Xcl and 
Xc2 production amplitudes. The previous assumption that the Tevatron events correspond 
to purely Xco CEP may therefore be unjustified. To give a prediction for the total cross 
section we must take into account the so-called 'enhanced' absorptive effects, which break 
the soft-hard factorization previously assumed in the evaluation of the Xco CEP cross 
section in [5] . The generalisation of the simple two channel eikonal model to include these 
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Figure 7: Xco, Xci and \ci differential cross sections, without experimental cuts, as a function of 
the photon p± and cos 9, where 9 is the angle between the [i + momentum in the J/ip rest frame 
and the direction of the Lorentz boost from the J/tp rest frame to the lab frame. Integrated cross 
sections normalised to unity before cuts are imposed in all cases. 



enhanced rescattering effects is outlined in [35], where the effect of including both eikonal 
and enhanced screening corrections, as well as non- forward outgoing protons in the hard 
matrix element, can be roughly accounted for by the introduction of an 'effective' survival 
factor Sgg. Recalling (1.1), it is found that the effect of enhanced absorption for the CEP 
of the light Xc is quite strong (see also [36, 51, 52]). In particular, ignoring for simplicity 
any impact parameter b dependence, we have 

(S 2 cS ) « (S 2 cnh ) x (S 2 cik ) « 1 (St ik ) , (5.8) 

with the enhanced survival factor S^ nh ( a t our present level of understanding) approxi- 
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Figure 8: Differential cross sections as a function of cos 0, with experimental cuts, and the 

p±, without experimental cuts. Integrated cross sections normalised to unity before cuts are imposed 

in all cases. 



mately the same for all three Xc states. Multiplying by this suppression factor and includ- 
ing the contribution from the three J states to the observed cross section, we can then use 
the branching ratio Br(x c o — > J/ip + l), as was done for the CDF data, to produce a rough 
value for the 'Xco' cross section at Tevatron energies (prior to any corrections due to the 
varying experimental acceptances of the Xc states), 

«65nb. (5.9) 

Vx=0 

Here we stress that a sizeable proportion of the observed events are predicted to correspond 
to Xd and Xc2 CEP. Thus, the combination of including the more general enhanced rescat- 



Xc 
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tering effects, which leads to a reduction in the predicted rate, and the contribution of the 
higher spin \ci and Xc2 states, which leads to an increase in the predicted rate, leaves the 
perturbative prediction for the total Xc cross section at the Tevatron largely unchanged 
and, crucially, still in good agreement with the experimental data. 

At this point we need also to consider the non-perturbative contribution to the cross 
section. Leaving the explicit evaluation for future work, we simply note that the previous 
calculations of [53, 54] show that the non-perturbative contributions of the three \ c states 
are also comparable, with the general results of Section 2 suggesting that, at least for 
the XcO and Xci, they will have a similar p± dependence to that outlined in Section 4. 
We can therefore reasonably assume that the relative values of (5.7) are approximately 
correct. However, given the overall uncertainty in the perturbative and non-perturbative 
cross section calculations, we note that the precise ratio of the Xc cross sections cannot be 
stated with certainty at this time. While our results suggest that some fraction of the 65 
± 10 candidate '^co' events observed at CDF are in fact Xci or Xc2 events, we can make 
no definitive prediction for their precise relative contributions. 

In fact, we have noted that the perturbative contribution to Xc2 CEP is strongly sup- 
pressed due to the J z = selection rule, for which it decouples from two real gluons. 
However, for the non-perturbative contribution we have no such selection rule, and we 
therefore cannot exclude the possibility that the Xc2 non-perturbative contribution is dom- 
inant. The above-mentioned Pomeron models suggest that this may be the case, although 
there remains a large degree of uncertainty in how to perform these calculations, and in 
particular which model of the Pomeron to choose. 

Putting aside the question of normalisation, we might hope to be able to distinguish 
between the three states by studying the angular and kinematical distributions of the final 
state particles as modelled in SuperCHIC, which should not depend strongly on the overall 
production rate. In Fig. 7 we show the polar angular distributions of the (i + p~ pair and 
the p± distributions of the photon for the three x states, as given by SuperCHIC. For the 
p± distribution we can see the clear separation in the x masses coming from the position 
of the Jacobian peaks as p± approaches the photon energy ~ M x — M^, although in 
the current experimental set-up we know it is not possible to resolve this separation. 

The angular distribution is more interesting: there is clearly a significant difference 
between the XcO and the Xci and Xc2 cases with, as expected from helicity conservation, 
the XcO decaying into purely transversely polarized J/Y>'s, while this is not the case for the 
Xci and Xc2- This in principle provides a way to determine if Xci and Xc2 mesons are being 
produced, irrespective of the particular mass resolution of the experiment. Unfortunately, 
this does not appear to be the case in practice, as we have yet to include the experimental 
cuts on the kinematics of the final state particles. At CDF we require in particular that 
the muon pseudorapidity \rj\ < 0.6 and the muon p± > 1.4 GeV. In Fig. 8 we show the 
angular distribution as before but with these cuts introduced, and immediately we can see 
that the clear difference in shape has not survived. We also show the p± distribution of 
the muons, from which it is clear that a sizeable fraction of the events will not pass the 
cuts (recalling that the cut must be passed by both muons), and the pseudorapidity cut 
further enhances this effect. Moreover, in the low (high) 8 region the p> + (~^ will be directed 
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along the motion of the J/ip with high p±, while the will be directed against the 

motion of the J/ip and will therefore have low p±. This is clear from Fig. 8 where the 
cos 6 w +1,-1 events which would have allowed us to distinguish between Xci,2 an d XcO 
production have not been accepted. We note that other potentially interesting variables, 
such as the difference in the azimuthal angle of the muons A</>^, appear to be equally- 
unpromising. It therefore seems that it will be very hard, given the experimental set up 
and the low statistics available at the present time, to distinguish between the three Xc 
states via the experimentally considered decay chain, although with more detailed analysis 
and/or higher statistics this conclusion may change. Another potential solution to this 
issue could be to consider a different decay chain, for example the direct decay of the Xc 
to charged hadrons, etc. (see Section 6 below). 

Finally, we can also see that the Xcl has a slightly higher acceptance, with in particular 
(at y x = 0) 

^^.cuts ^^.cuts ^^.cuts 

-TTT ■ -TTT ■ -TTT ~ 16 % : 17 % : 15 % ■ ( 5 - 10 ) 
dcr tot dcr tot d(T tot 

The acceptance is therefore reasonably uniform, and so should not present a significant 
obstacle when considering Xcl and Xc2 production experimentally. 



6. Summary and Outlook 

Motivated by the recent experimental observation of exclusive Xc events at the Tevatron, we 
have updated the earlier studies of central exclusive scalar XcO meson production to include 
Xcl and Xc2 mesons. Due to the low scale, M Xc /2, and very large rapidity gap coverage 
(At? ~ 7.4 units) in the CDF measurement [20], the contamination from processes in which 
the incoming protons dissociate is relatively small. The CDF Xc event selection therefore 
effectively ensures that they come from the exclusive reaction, pp — > p + Xc + P- 

Although XcO production was previously assumed to be dominant, we find that the XcO, 
Xcl and Xc2 rates for the experimentally considered Xc — ► J 1^1 — * M + / u_ 7 decay process 
are in fact comparable. We have developed a new Monte Carlo event generator, Super- 
CHIC, which models the central exclusive production of the three Xc states via this decay 
chain, and have used this to explore possible ways of distinguishing them, given that their 
mass differences are not resolvable within the current experimental set-up. Although we 
find that the severity of current experimental cuts appears to preclude this discrimination, 
the acceptance does not change crucially between the three states and so our conclusions 
regarding the overall rates remain unchanged. This therefore raises the interesting possi- 
bility that exclusive Xcl and Xc2 production has already been observed at the Tevatron. 
Higher statistics and/or a broader acceptance coverage for the photon and leptons could 
help discriminate between the Xc mesons via differences in the angular correlations between 
the final-state particles. We note also that the addition of forward proton detectors would 
certainly allow discrimination between the different C-even states via the measurement of 
the relative azimuthal angular distribution between the outgoing protons. 

To further resolve the spin-parity assignment issue in the absence of forward proton 
detectors, it would be instructive to observe central exclusive Xc production in other decay 
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channels, in particular nir or KK, see [5]. These modes are ideally suited for spin-parity 
analysis: the -kit or KK decay modes of the XcO meson have a branching fraction of about 
1%, while these decay channels are forbidden for Xd and suppressed by about a factor 
of 5 for the Xc2 relative to the XcO> i n contrast to the \c2 relative enhancement for the 
J/tpj channel. Another interesting mode for discriminating between the CEP of different 
Xc states is Xc — > PP, since the branching fraction for XcO (— 0.024%) is a factor of 3 higher 
than that for xa,2 [37]. The AA mode (branching fraction for XcO — 0.034%) could also be 
important for spin-parity analyzing. 18 

In the case of two-body final states (tttt, KK, pp) where the very forward protons 
are not detected, further cuts can be imposed to reduce the contribution of events where 
the protons dissociate (via single and double diffractive dissociation). These include, for 
example, cuts on the transverse momentum of the resonance and on the final particles' 
accoplanarity angle (in the frame where the rapidity of the resonance is zero). 19 

In this paper we have focused on Xc meson production at the Tevatron. It is of course 
straightforward to extend our results to the LHC, and we will consider this is in a future 
study [57]. Note that we do not expect the Xc CEP rate to have a strong energy dependence 
when going from the Tevatron to the LHC. The growth of the bare amplitude caused by 
the increase in the gluon density at smaller x is compensated by a smaller gap survival 
factor at the larger LHC energies, especially S 2 nh , the value of which decreases due to the 
larger rapidity interval available for the 'enhanced' absorptive corrections, see Fig. 1(d). 
Indeed, the measurement of the ratios of the CEP rates at the two different (Tevatron 
and LHC) collider energies could allow the effects of enhanced absorption to be probed, 
since in these cross section ratios various uncertainties (for example, NLO corrections to 
the gg — > x transition etc.) would cancel out. 

The issue of forward proton detection is more relevant at the LHC (in particular 
for Xb states) as the planned near-beam proton detectors, see Refs. [11, 15, 58], would 
allow us to measure the outgoing very forward protons. As we have already noted, the 
azimuthal angle distributions would provide interesting additional information with which 
we could discriminate between different J PC states as well as investigate the dynamics of the 
survival factors S 2 . Moreover, as pointed out in [40], measuring the transverse momentum 
and azimuthal angle correlations between the outgoing protons would allow us to probe 
the proton opacity Q(s,bt) and perform a detailed test of the whole diffractive formalism. 
Even before the forward proton detectors become operational, a broad programme of heavy 
quarkonium studies can be performed with the existing LHC detectors (ALICE, ATLAS, 
CMS and LHCb), especially if the rapidity gap coverage is increased by using forward 
shower counters (FSC) along the beam line [59], thereby allowing detection and triggering 
on rapidity gaps in diffractive events. 

In future work we will also consider Xb and r]b,c meson CEP and include them in the 

18 Rough estimates show that the background from continuum it-k, KK, pp and AA central production 
should be quite manageable. This is in accord with measurements in two-photon collisions [55], where the 
XcO, 2 resonances decaying to tt-k and KK final states are clearly seen. 

19 This procedure is similar to that used in the separation of exclusive lepton-pair production via photon- 
photon fusion, pp — > p + l + l~ + p, see [56] for more details. 
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Monte Carlo. The perturbative Xbo CEP rate was calculated in [3, 5] analogously to the 
XcO- While the higher M x scale results in a more perturbatively reliable prediction, the 
reduced rate suggests that exclusive Xb production, observed through the Xb — ► T(1S)7 — > 
/i + /i~7 decay chain, may only be relevant at the LHC. The experimental values for the 
Xbj — ► T(1S)7 branching ratios carry larger uncertainties, but a similar hierarchy to the 
Xc case is observed, and so the contribution of higher spin states should once again be 
considered, with in particular the result of (3.20) for the relative rates remaining valid. 
However, we note that the higher mass scale will result in a stronger suppression of the 
higher spin states relative to the Xbo than in the Xc case. In particular, for Xbi production 
we have an explicit factor of in the denominator of (3.20), while the larger expected 
value of (Qj_) will result in a stronger suppression of the Xb2 state. 

The rj Cj i, CEP cross sections can be calculated by using the same formalism as for the 
X mesons, the only difference being that for the L = quarkonium state the gg — > r\ vertex 
is proportional to the value of the wave function (f)(0) at the origin and not to 4>'(0) as in 
(3.12). The vertex should therefore be normalised to the leptonic width of the J/ip (or 
T for rjb) decay - the vector mesons from the same L = multiplet. Note also that for the 
heavier Xb or r/b mesons we expect a slightly lower slope b eS (due to a larger mean (Q\)) 
and a larger value of S^ nh (due to a smaller rapidity interval available for the 'enhanced' 
absorptive corrections, Fig. 1(d)). Preliminary estimates indicate that the r/ c CEP rate 
(which we recall from [4] will be proportional to {p( p% } / (Qj_) 2 ) ls expected to be about 
two orders of magnitude lower than in the XcO case. 

We are also planning to revisit 77 CEP in a wider interval of photon Et, rapidity 
and di-photon mass M than that considered in [23] and to include this in the Monte Carlo 
generator. Note that the measurement of the ratio of 77 CEP at Et = 5 GeV to that 
of Xb production may allow us to reduce various uncertainties in the calculations, with in 
particular the dependence on the survival factors cancelling out. 

Finally, we note that the spin-parity analyzing properties of central exclusive produc- 
tion could shed light on the dynamics of the zoology of 'exotic' charmonium-like states 
(X,Y,Z) which have been discovered in the last few years (see for example [37, 60]), and 
whose nature and in many cases spin-parity assignment still remain unclear. 
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A. Xc — > <?<? amplitudes 
A.l Xco 

We use the formalism of [44] and the kinematics of Section 3 throughout. The general 
colour-averaged vertex for the coupling to two gluons has the form 

V ee ^4vSp = \Jl^i(M 2 + (gift)) " 2/ 2 °) , (A.l) 
where ei )2 are the polarization vectors of the incoming gluons and 

qi = xipi + Q± - Pi ± , (A.2) 

<?2 = X 2 P2 ~Q±- P2 ± , (A.3) 



while we define 



I°i = F^F 2 ^ , (A.4) 
7° = qlF^F 2 >^q 2 , a . (A.5) 



Here F^ u is the usual field strength tensor for the gluons and 

1 4 5s 2 re 



2VNc (QiQ2) 2 V 4vrM 



0c(O) , (A.f 



where g s is the strong coupling and 4>' c (0) is the derivative of the Xc radial wavefunction at 
the origin. Note that our definition of c differs from that of Ref. [44] by a factor 

xab 1 

(M;3k\l)tU b ik = — -= -► — = , (A.7) 

where (3i; 3A;|l) is the colour space Clebsch-Gordon coefficient for the colour singlet quark 
configuration, and we have averaged over the gluon colour indices a, b in the last step. 
Recalling (3.2), we make the replacement e^e^V^ — > -pIp^V^ and use the gauge invariance 
of Vpv to give 

A° = 2(gi ±(?2± ) , (A.8) 
I°2 = qUL • (A.9) 

Here we have made use of the identity 

(qm) = \{M 2 - q 2 ± - q 2 J . (A.10) 



We therefore obtain 



Vo = \flw {{qi ^ ){3M * - ^ " - 2 ^L) ( A - n ) 
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A.2 Xd 

The general colour-averaged vertex has the form 



Vi = --(ll+lh, (A.12) 



where 



I\ = e^efF^qJ (A.13) 
II = e^efF^qJ . (A.14) 

Here is the xi polarization vector and ^ VOL ^ is the antisymmetric Levi-Civita tensor. 
This gives 

lie 

Vi = -—e^ef Pl , uP2 , a ((q 2± Uq 1± ) 2 - (q 1± ),(q 2± f) • (A.15) 

A.3 X c2 

The general colour-averaged vertex has the form 

V 2 = -cV2MI% , (A. 16) 

where 

if = zT F f F l[> ■ (A-17) 
e^ a is the xi polarization tensor, which satisfies 



— ^v^i ■, tjj 1 — , e^ u P£ — , (A. 18) 



Yl 6 ^ e «/3 = |( + V^Pua) ~ I V P "/3 , (A-19) 

pol 

pi 1 pv 

W v _ _ gt xv + _xx _ ( A _20) 



We therefore obtain 



V 2 = ^^(s(qn.)^q 2 ±) a + 2(q 1± q 2± ) Pl ^ 2a )e* x ^ . (A.21) 



B. Xc and J/t/> decay amplitudes 

We will make use of the identities 



= -4##, (B.l) 
e^ a/3 W = -3!4 M ^^ > (B-2) 
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B.l Xc(0 ++ ) - J/il> + 7 



The scalar xo decays into a transversely polarized photon with a uniform angular distribu- 
tion in its rest frame, and conservation of angular momentum therefore requires the J/ip 
to be transversely polarized. 



B.2 Xc (l ++ ) - J/1> + 7 

The amplitude that is expected to dominate (as it is the amplitude which corresponds to 
the dipole transition [44]) is of the form 

M ~ ^eXefple}* . (B.4) 

Squaring and summing over photon polarizations, and making use of (B.2), we find 

^|^ 1 | 2 ~|(p 7 e v ,)| 2 + |(p 7 e x )| 2 + 2Re[(e x e^)(e*p 7 )(p 7 e;)] . (B.5) 

€-, 

The normalisation is given by summing over J/tp polarizations and making use of (B.l) 
and (B.3) 

\Ar m \ 2 -\(e xP ,)\ 2 + ^^(( W ) + 2Re[(6> 7 )(6 x ^)]) . (B.6) 

We then divide by the normalisation factor to give the relative amplitudes squared for the 
three different xi polarizations. 



B.3 Xc(2 ++ ) - J/ip + 7 

Following similar arguments to the xi case (that is, assuming the dipole transition domi- 
nates), we can write the invariant amplitude as 

A 2 = e^^f^U ■ (B.7) 

To study angular correlations it is necessary to consider the explicit form of the X2 polariza- 
tion tensor. This represents the irreducible tensor operator for J = 2 angular momentum, 
which can be decomposed in terms of the spin and orbital polarization vectors [44] 

eil z) = ^2ei Sz) e^(S = l,L = l,S z ,m\J = 2,J z ) , (B.8) 

S z ,m 

where (S = 1, L = 1, S z ,m\ J = 2, J z ) are the Clebsch-Gordon coefficients and e^ z \ th^ 
have the usual explicit representation in (say) the X2 rest frame. We can thus decompose 
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the 5 polarization states as 




(B.9) 



(V+ e o i e o e +) 

(e+e- +2e°e°+e-e+) 
f e - e o + e o e -~) 



(B.10) 



(B.11) 



(B.12) 



(B.13) 



Returning to (B.7), we obtain 



- eV{g uu ((PyP^t - (p^)pi)((p lPi ,)et - (p^pf) 
+ pl P l(M*e* c ?et-pipt)), 



(B.14) 



where the normalisation is given by 



l^ om f ~ V" (Gw) V + M^pjjpl) + 2( W )pX Re[^ Q e* a " ] • (B.15) 



References 

[1] D. Robson, Nucl. Phys. B130 (1977) 328; 
F.E. Close, Rept. Prog. Phys. 51 (1988) 833. 

[2] P. Minkowski, Fizika B 14 (2005) 79 [arXiv:hcp-ph/0405032]. 

[3] V. A. Khoze, A. D. Martin and M. G. Ryskin, Eur. Phys. J. C 19, 477 (2001) [Erratum-ibid. 
C 20, 599 (2001)] [arXiv:hep-ph/0011393]. 

[4] A. B. Kaidalov, V. A. Khozc, A. D. Martin and M. G. Ryskin, Eur. Phys. J. C 31, 387 
(2003) [arXiv:hep-ph/0307064]. 

[5] V. A. Khozc, A. D. Martin, M. G. Ryskin and W. J. Stirling, Eur. Phys. J. C 35, 211 (2004) 
[arXiv:hep-ph/0403218]. 

[6] A. Kaidalov et al, V.A. Khozc, A.D. Martin and M. Ryskin, Eur. Phys. J. C 33 (2004) 261, 
hep-ph/0311023. 

[7] E. Klcmpt and A. Zaitsev, Phys. Rcpt. 454 (2007) 1 [arXiv:0708.4016 [hcp-ph]]. 

[8] F. E. Close and A. Kirk, Phys. Lett. B 397 (1997) 333 [arXiv:hcp-ph/9701222]; 

F. E. Close, A. Kirk and G. Schulcr, Phys. Lett. B 477 (2000) 13 [arXiv:hcp-ph/0001158]. 

[9] S. Heinemeyer, V. A. Khoze, M. G. Ryskin, W. J. Stirling, M. Tasevsky and G. Weiglein, 
Eur. Phys. J. C 53 (2008) 231 [arXiv:0708.3052 [hcp-ph]]. 

[10] V. A. Khozc, A. D. Martin and M. G. Ryskin, Eur. Phys. J. C 14, 525 (2000) 
[arXiv:hcp-ph/0002072]. 



-27- 



[11] M. G. Albrow and A. Rostovtsev, arXiv:hep-ph/0009336. 

[12] V. A. Khozc, A. D. Martin and M. G. Ryskin, Eur. Phys. J. C 23, 311 (2002) 
[arXiv:hep-ph/0111078]. 

[13] A. De Rocck, V. A. Khoze, A. D. Martin, R. Orava and M. G. Ryskin, Eur. Phys. J. C 25, 
391 (2002) [arXiv:hcp-ph/0207042]. 

[14] For a recent review see A. D. Martin, M. G. Ryskin and V. A. Khoze, arXiv:0903.2980 
[hep-ph]. 

[15] M. G. Albrow et al. [FP420 R&D Collaboration], arXiv:0806.0302 [hep-ex]. 

[16] P. Bussey and P. Van Mechelen in: H. Jung et al, arXiv:0903.3861 [hep-ph], p. 557. 

[17] C. Royon, Acta Phys. Polon. B 39, 2339 (2008) [arXiv:0805.0261 [hep-ph]]. 

[18] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 99 (2007) 242002 [arXiv:0707.2374 
[hep-ex]]. 

[19] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. D 77,(2008) 052004 [arXiv:0712.0604 
[hep-ex]]. 

[20] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 102, 242001 (2009) 
[arXiv:0902.1271 [hep-ex]]. 

[21] M.Albrow, arXiv:0909.3471 

[22] Mike Albrow and Jim Pinfold, private communication. 

[23] V. A. Khoze, A. D. Martin, M. G. Ryskin and W. J. Stirling, Eur. Phys. J. C 38 (2005) 475 
[arXiv:hcp-ph/0409037]. 

[24] J. Pumplin, Phys. Rev. D47 (1993) 4820. 

[25] F. Yuan, Phys. Lett. B 510, 155 (2001) [arXiv:hep-ph/0103213]. 

[26] V. A. Petrov and R. A. Ryutin, JHEP 0408 (2004) 013 [arXiv:hcp-ph/0403189]; 
V. A. Petrov, R. A. Ryutin, A. E. Sobol and J. P. Guillaud, JHEP 0506 (2005) 007 
[arXiv:hep-ph/0409118]. 

[27] A. Bzdak, Phys. Lett. B 619 (2005) 288 [arXiv:hcp-ph/0506101]. 

[28] M. Rangel, C. Royon, G. Alves, J. Barrcto and R. B. Peschanski, Nucl. Phys. B 774, 53 
(2007) [arXiv:hep-ph/0612297]. 

[29] R. S. Pascchnik, A. Szczurck and O. V. Teryaev, Phys. Lett. B 680 (2009) 62 
[arXiv:0901.4187 [hep-ph]] and references therein. 

[30] G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. D 51 (1995) 1125 [Erratum-ibid. D 
55 (1997) 5853] [arXiv:hep-ph/9407339]. 

[31] N. Brambilla, A. Pineda, J. Soto and A. Vairo, Rev. Mod. Phys. 77 (2005) 1423 
[arXiv:hep-ph/0410047]; 

N. Brambilla and A. Vairo, Acta Phys. Polon. B 38, 3429 (2007) [arXiv:0711.1328 [hep-ph]]. 

[32] V. A. Khoze, A. D. Martin and M. G. Ryskin, Eur. Phys. J. C 18, 167 (2000) 
[arXiv:hcp-ph/0007359]. 



- 28 - 



[33] M. G. Ryskin, A. D. Martin and V. A. Khoze, Eur. Phys. J. C 54, 199 (2008) 
[arXiv:0710.2494 [hcp-ph]]; 

M. G. Ryskin, A. D. Martin and V. A. Khozc, Eur. Phys. J. C 60 (2009) 249 
[arXiv:0812.2407 [hcp-ph]]. 

[34] V.A. Khozc, A.D. Martin and M.G. Ryskin, JHEP 0605, 036 (2006) [arXiv:hcp-ph/0602247]. 

[35] M. G. Ryskin, A. D. Martin and V. A. Khoze, Eur. Phys. J. C 60 (2009) 265 
[arXiv:0812.2413 [hcp-ph]]. 

[36] J. Bartels, S. Bondarenko, K. Kutak and L. Motyka, Phys. Rev. D 73 (2006) 093004 
[arXiv:hcp-ph/0601128]. 

[37] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008) and 2009 partial update 
for the 2010 edition. 

[38] L.D. Landau, Dokl. Akad. Nauk SSSR 60 (1948) 213; 
C.N. Yang, Phys. Rev. 77 (1950) 242. 

[39] D. Barberis et al. [WA102 Collaboration], Phys. lett. B440 (1998) 225; ibid. B422 (1998) 
399; 

A. Kirk et al. [WA102 Collaboration], arXiv:hcp-ph/9810221. 

[40] V. A. Khoze, A. D. Martin and M. G. Ryskin, Eur. Phys. J. C 24, 581 (2002) 
[arXiv:hcp-ph/0203122]. 

[41] S. Aid et al. [HI Collaboration], Nucl. Phys. B 472 (1996) 3 [arXiv:hep-ex/9603005]. 

[42] A. D. Martin and M. G. Ryskin, Phys. Rev. D 64, 094017 (2001) [arXiv:hcp-ph/0107149]. 

[43] A. G. Shuvaev, K. J. Golec-Biernat, A. D. Martin and M. G. Ryskin, Phys. Rev. D 60, 
014015 (1999) [arXiv:hcp-ph/9902410]. 

[44] J. H. Kuhn, J. Kaplan and E. G. O. Safiani, Nucl. Phys. B 157, 125 (1979). 

[45] R. Barbicri, R. Gatto and R. Kogcrlcr, Bound Phys. Lett. B 60, 183 (1976). 

[46] V. A. Novikov, L. B. Okun, M. A. Shifman, A. I. Vainshtein, M. B. Voloshin and 
V. I. Zakharov, Phys. Rcpt. 41, 1 (1978). 

[47] F. E. Close, G. R. Farrar and Z. p. Li, Phys. Rev. D 55, 5749 (1997) [arXiv:hcp-ph/9610280]. 

[48] R. Barbicri, M. Caffo, R. Gatto and E. Rcmiddi, Phys. Lett. B 95, 93 (1980). 

[49] R. Alkofer and C. S. Fischer, Fizika B 13, 65 (2004) [arXiv:hep-ph/0309089]. 

[50] M. Gluck, E. Rcya and A. Vogt, Z. Phys. C 67, 433 (1995). 

[51] L. Frankfurt, C.E. Hyde, M. Strikman and C. Weiss, Phys. Rev. D75, 054009 (2007); 
arXiv:0710.2942 [hcp-ph]; 

M. Strikman and C. Weiss, arXiv:0812.1053 [hep-ph]. 

[52] E. Gotsman, E. Levin, U. Maor and J. S. Miller, Eur. Phys. J. C 57, 689 (2008) 
[arXiv:0805.2799 [hcp-ph]]. 

[53] H. A. Peng, Z. M. He and C. S. Ju, Phys. Lett. B 351, 349 (1995). 

[54] E. Stein and A. Schafer, Phys. Lett. B 300, 400 (1993). 



-29- 



[55] T. Mori et al. [BELLE Collaboration], J. Phys. Soc. Jap. 76 (2007) 074102 [arXiv:0704.3538 
[hep-ex]]; 

H. Nakazawa et al. [BELLE Collaboration], Phys. Lett. B 615 (2005) 39 
[arXiv:hep-ex/0412058] ; 

S. Uehara et al. [BELLE Collaboration], arXiv:0903.3697 [hep-ex]. 

[56] V. A. Khozc, A. D. Martin, R. Orava and M. G. Ryskin, Eur. Phys. J. C 19, 313 (2001) 
[arXiv:hcp-ph/0010163]. 

[57] L. A. Harland-Lang et al, in preparation. 

[58] M. Albrow et al, CERN-LHCC-2006-039, CERN-LHCC-G-124, 
CERN-CMS-NOTE-2007-002, Dec 2006; 

G. Anclli et al. [TOTEM Collaboration], JINST 3 (2008) S08007. 

[59] M. Albrow et al. [USCMS Collaboration], arXiv:081 1.0120 [hep-ex]; 
J. W. Lamsa and R. Orava, arXiv:0907.3847 [physics. acc-ph]. 

[60] G. V. Pakhlova, arXiv:0810.4114 [hep-ex]; 

Phys. Atom. Nucl. 72 (2009) 482 [Yad. Fiz. 72 (2009) 518]. 



- 30 - 



